Inverse quantum Zeno effect and power laws in a model chiral molecule 
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We consider a chiral molecule modeled by two sets of rotational internal energy levels of different 
parity and degenerate ground states coupled by a constant interaction. By assuming that the 
collisions in a gas of identical molecules, provoke transitions between adjacent energy levels of 
the same parity, we apply the prescriptions of the continuous time quantum random walk to a 
single molecule, interpreted as an open quantum system, and build up the master equation driving 
its internal dynamics in case of a general distributions of collision times. The asymptotic exact 
dynamics of the coherence terms and the populations of the energy levels are analytically evaluated 
for relevant classes of non-Poissonian distributions of the collision times. The inverse quantum Zeno 
effect and inverse power law relaxations emerge over estimated long time scales. 
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I. INTRODUCTION 

The Quantum Zeno effect (QZE), originally introduced 
by Misra and Sudarshan [ij in the early 1977, is gener- 
ally interpreted as the hindrance of the dynamics of an 
unstable quantum system, caused by frequent measure- 
ments. Experimental evidence of the QZE were shown 
by Cook [2] in the early 1988, by Itano et al. in forced 
Rabi oscillations between discrete atomic levels, and in 
spontaneous decaying systems by Fisher et al. [3], to 
name a few. Under particular conditions the QZE is also 
obtained over short time scales 0-01, the literature on 
this argument is wide. The inverse quantum Zeno effect 
(IZE) is, on the contrary, the acceleration of the relax- 
ation process due to repeated measurements. This effect 
has been widely studied through theoretical approaches 
d, Q and experimental evidences Q- In Ref. @ the 
transition from the QZE to the IZE is obtained by chang- 
ing the time interval between successive measurements. 
Also, inverse power law behavior emerge over long time 
scales in the time evolution of the survival probability of 
an unstable quantum system [lo| . 

Recently [ll|, a gas of colliding identical chiral model 
molecules has been considered and the Rabi-like oscil- 
lations between the two sets of rotational levels of the 
model molecule are shown to be inhibited by an increase 
of the mean collision time. This behavior has been inter- 
preted as the appearance of the QZE. This picture refers 
to the experimental behavior of the nuclear spin depolar- 
ization in ^^CH^F molecule [l^ . The quantum number 
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of the total spin of the three protons takes the values 3/2 
(orto) and 1/2 (para), the transition between states of 
different parity are forbidden in the electric dipole inter- 
action and the spin flip emerges from a weak coupling 
between two levels of different spin parity. Here, the 
freezing of the spin relaxation due to an increase of the 
gas pressure, is considered as the appearance of the de- 
phasing caused by molecular collisions. 

In this scenario, we aim to study the dynamics of a 
chiral molecule modeled in Ref. [il| by considering rel- 
evant cases of non-Poissonian distributions of collision 
times and recover the Poisson statistics as a particular 
case. The construction of the master equation is per- 
formed through the prescription of the continuous time 
random walk (CTRW), originally developed by MontroU 
and Weiss [1^. If the distribution of time intervals be- 
tween two consecutive jumps is non-Poissonian, anoma- 
lous diffusion is obtained from the random walk. This 
technique is adapted to the dynamics of an open quan- 
tum system [lj| randomly interacting with its external 
environment. In case the distribution of the interval be- 
tween two consecutive interaction is non-Possonian, the 
master equation mimicking the dynamics of the system 
of interest results to be non-Markovian. This procedure 
is also named as the continuous time quantum random 
walk (CTQRW) and is developed in Refs. ilB] and [ll. 
This powerful technique recovers a non-Markovian mas- 
ter equation as a subordination [l3| to a Lindblad dy- 
namics [H, [l^. Once the master equation is built up, 
we study the exact time evolution and the eventual ap- 
pearance of the IZE over long time scales. 

Details on the construction of the master equation cor- 
responding to a general distribution of collision times are 
given in Section [TTl Section |IlT] is devoted the study of 
the asymptotic dynamics of the populations of the energy 



2 



levels for a general distribution of collision times and rel- 
evant particular cases. Time scales for inverse power laws 
are analytically estimated and the appearance of the IZE 
is discussed. In Appendixes |X] and [B] a detailed analysis 
of the convoluted structure equations driving the dynam- 
ics of the levels populations and the decoherence terms 
is performed. 



II. THE MASTER EQUATION FOR 
NON-POISSONIAN DISTRIBUTIONS OF 
COLLISION TIMES 

The model multilevel molecule is considered as an open 
quantum system and the collisions with other molecules, 
represented as an interaction with an external environ- 
ment, are described by the action of a set of superoper- 
ators on the density matrix of the molecule itself. This 
picture will be adopted for the construction of the master 
equation mimicking the dynamics of the internal energy 
levels of the model molecule in case of a general distri- 
bution of the collision times. By following the procedure 
of the CTQRW, a non-Markovian master equation is ob- 
tained for the time evolution of the internal energy levels 
in case of non-Poissonian distributions of collision times. 
Recent applications of this technique are shown in Ref. 
[20] in order to interpret the fluorescence fluctuations in 
blinking quantum dots. 

Briefly, the CTQRW provides the following time evo- 
lution of the density matrix: 



Pit) 



n=0 



Voit-t') W„(t')[p(0)] dt'. (1) 



In the CTRW, the function Wnit) recalls the probabil- 
ity that n jumps has occurred, the last one at time t, 
while the function Vo{t) recalls the probability that no 
jumps occurs in time interval t. For our model, Vo{t) and 
Wn (t) are superoperators properly defined through their 
Laplace transforms [2J]: 

Voiu)^ Poiu^Co), Wniu) ^ (Aw {u - Co)r , (2) 

where the function Po{t) is the probability that no colli- 
sion occurs in a time interval t, while the function w{t) 
is the statistical distribution of the random intervals be- 
tween two consecutive collisions. The superoperator A, 
mimicking the effect of each collision, is related to the 
Liouvillian Ci by the relation: Cj — A — I, where / 
is the identity superoperator. In this way, the Laplace 
transform of Eq. ([T|) gives: 



w{u- Co) 



u - 







Aw {u — Co) 



[pm 



(3) 



equivalent to the following non-Markovian master equa- 
tion: 



p{t) ^ Co W)] 



^{t-t')Ci 







M'-'') [p{t')^ dt', 

(4) 



where $ = uw{u)/ (1 — w{u)) is the Laplace transform of 
the memory kernel related to the distribution of collision 
times. For more details we refer to [l^ HgI. [2l|. 

We now remind the model adopted in Ref. 11]: the 
spatial dynamics is neglected and each molecule is mod- 
eled by the internal rotational energy levels of different 
parity, labeled as left (L) and right {R). The ground 
levels are energy degenerate and coupled by a constant 
interaction, the collisions cause transitions between ad- 
jacent energy levels of the same parity, while the transi- 
tions between L and R spin states are forbidden, except 
through the ground states. The model also provides that 
the molecule does not reach the dissociation energy. For 
more details, we refer to 11]. 

Briefly, the Schrodinger equation in Ito form, corre- 
sponding to a Poisson distribution of collision times of 
mean tq, gives the following master equation describing 
the time evolution of the statistical density matrix in the 
Hilbert space of the total spin: 



m = [Cf + ^^c.] [pit)] 



(5) 



The Liouville superoperator L f drives the free evolution, 
Cf [■] — —i/h[H, ■] and the total Hamiltonian is deflned 
as follows: 



H — Ho + Hi, 



(6) 



E 

n_L = l 



Nr 



Ho^ y ] Eni^\nL){nL\ + ^ EnJnR){nR\, 

nR = l 

Hi=hni\lL){lR\ + \lL){lB\), 

Ei,=Ei„=Ei, (U||1;^)=0. 

The term Hi mimics the constant interaction between 
the ground states of different parity, while the superop- 
erator Cc, related to the action of each collision, reads as 
follows: 



[.] = ^^[V,■]~^[V[V,■]], 



(7) 



a.s {\ns){ns+i\ + \ns + l){ns\) ■ 



s=L,R ns = l 

The structure of the interaction Hamiltonian V recovers 
the assumption that the collisions provoke transitions be- 
tween "nearest neighbor" energy levels of the same parity. 
In the Poisson case, the distribution of collision times and 
the corresponding memory kernel, respectively, read 

. ^ e-*/^° ^ , , 5U) 
wpit) = , $p t - — 

To To 

and the Lindblad structure [3, [H of Eq. © IS recov- 
ered from Eq. Q. The master equation ^ describes 
the dynamics of the multilevel molecule undergoing the 
collision interactions © with a generic distribution of 
the collision times, provided that Co = Cf and Cj = Cc- 
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This way, the following non-Markovian master equation 
is obtained: 

p{t) = C}p{t) + f ^{t- t') [e^4*-*') [p {t')]\ dt'. 
Jo 

(8) 

Let To be the characteristic time of the Poisson distri- 
bution of collision times and AE be the smallest energy 
difference between every couple of L and R states, the 
ground states are excluded. The following estimates: 
AE/h ~ 10~^s, ~ IkHz and tq ~ 1/is, justify the 
assumption that all the level pairs are far from resonance 
condition. We also assume that the investigated physical 
conditions fulfill the constraint 



AE 1 
— — ^ max <\l, — 



(9) 



holding true for every characteristic time t$ of each dis- 
tributions of collision times that we are going to consider, 
being the corresponding memory kernel. As pointed 
out in Ref. the inequality (jH]), reveals the existence 
of a "fast" dynamics giving a vanishing average contribu- 
tion over long time scales, this issue is discussed in details 
in Appendix 1^1 where the master equation ([SJ is studied. 



III. TIME SCALES FOR THE INVERSE 
QUANTUM ZENO EFFECT AND POWER LAWS 

This Section is devoted to the dynamics of the popula- 
tions of the energy levels and to decoherence process. We 
aim to show the occurrence of the IZE and inverse power 
law relaxations over estimated long timescale, for rele- 
vant classes of non- Poisson distribution of collision times, 
recovering the Poisson statistic as a particular case. 

The coherence term, 

and the populations of the whole L and R levels. 



N 

E 



are evaluated in Appendix for a generic distribution of 
collisi on times, through Eqs. dM]), (HH), (pi2t . ((M3|) . 
(|A14[) and those obtained by exchanging the indexes L 
and i?, holding true for times t ^ 

The detailed study performed in Appendix [Bl reveals 
that the populations of the whole L and R energy levels 
are driven uniquely by the coherence term p'^it) and the 
populations of the ground levels: 

PUt) = np-{t), Pnit) = -nf{t), (10) 
p-{r)^2n{p,^{t)-p,,{t)). (11) 



The final value theorem [221, applied to Eqs. ([BT]), p32]) . 
([TU]) and pT|) . suggests the existence of stable asymptotic 
configuration for every distribution of collision times. 



Pl,,, (+oo) 



Pl (+oo) _ ol 
Pr (+oo) ur 



(12) 



The whole population of both the L and R levels tend 
to a stationary distribution depending uniquely on the 
collision interaction, represented by the parameters aL 
and Qfij, and is independent of either the Rabi-like oscil- 
lations between the ground levels or the statistical distri- 
bution of the random collision times. In case ul = ctR, 
the asymptotic populations of the whole L and R energy 
levels become identical. 

We now show in detail the exact dynamics of the lev- 
els populations and the coherence terms over long time 
scales, for relevant non-Poissonian distributions of colli- 
sion times. 



A. Fractional diffusion 

As a first case, we consider the fractional diffusion pro- 
cesses widely studied in literature [H, [H [H, Hi] . The 
corresponding distribution of collision times. 



Writ) = alt-^^ Ei^2rA-2r (-a^'"'") 



> r > 0, 



(13) 

is defined through the generalized Mittag-Leffler function 

MM 



z,/3eC, 5R{a}>0. 



We notice that the mean time is infinite and the Poisson 
statistic is recovered for r = 0. 

We plug the Laplace transform of the corresponding 
memory kernel. 



UWrju) 2 2r 

$r(M) = -; = a^u , 

1 — Wr(U) 



l>r>0, (14) 



in the general expression (jBl[) and analyze the behavior 
for It 0+. The Laplace inversion gives the following 
asymptotic inverse power law behavior: 



{aR - aL) t-^' 



2Vlar{aL +afl,) r(r- 1/2) 



t -> -l-oo. 



(15) 

Initially, only the ground L level is populated, which 
means Pl{^) = 1- Over a long time scale, t ^ Tr, defined 
below, Eqs. (|10p and (|Bip lead to an inverse power law 
relaxation described by the following asymptotic form: 



+ (-) 



aL + aR 
(aR - aL) t-^'^ 
2ar{aL+aR)'^T{r + 1/2)' 



-oo. (16) 



A possible choice of the long time scale arises from the 
constraints for the convergence of the series expansion of 
Eq. (jBip . In this way, the following form in dimension- 
less units is obtained: 
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X (5 + 2017^ + arQ;fl,(ll + 4417^ + 40^07^(8 + 12VL^ + a,.a7j)))^ + (l + 4f^^) (2 + QraR^'i + arai?X9 + Sa^aj?^ 
X (5 + 2arafl) ))) + 2a,. (l + 4f]2) aL(^2 + a^ai? (4 + arai?(9 + 20^0;^ (5 + 2a,.ai^) ))^ + (l + 417^) q;| 

X (^9 + 2a,.afl(9 + 2a,.a7?(l0 + a^aR (11 + 4a,.aj^) ))))/( 2^aln'^ala]i {ul + a^) 



2/(l-2r). 



(17) 



For large collision time intervals, t ^ ar , the 

corresponding statistical distribution is described by the 
asymptotic form 



, , (1 - 2r)i-2(i-r) 



-00. 



Thus, the following behavior emerges: a decreasing of 
the probability that two consecutive collisions occur at 
large time intervals, i.e., an increase of the parameters 
ttr, provokes, over long time scales, t ^ r^, a faster re- 
laxation of the initially populated L levels. The inverse 
power law relaxations, occurring over long time scales, 
are accelerated by a decreasing of the probability that 
two consecutive collisions occur at large time intervals. 
This correspondence could be identified as an IZE over 
long time scales, occurring for large collision time inter- 
vals. 



B. Inverse Power law distribution 

We now consider a power law distribution of collision 
times, 



(A^-i)r^' 

{t + Tf 



2 > /i > 1, 



(18) 



of infinite mean time. This case has been widely studied 
in literature in both classical [27| and quantum processes 



[20|. The Laplace transform of the distribution of colli- 
sion times 



w^{u)^l-T{2- ^i){uTY ^ + 



uT 
2-^' 



u -> 0+, 
(19) 

gives the following asymptotic form of the memory ker- 
nel: 



,2-M 



^'^^^^ - T^tR ^ Ti2-,)T.-^ - " - 

(20) 

The analysis of Eq. (IBip for m — 0+, gives inverse power 
law relaxations over a long time scale, t ^ t^, defined 
below, 



t 

T 



[aR - aL) yr (2 - ^) 
2nT{aL+aR)''T{{l- /2) 

-(1+m)/2 



t — > -f 00, 



(21) 



OIL + OLR 
^ ^^(l-p)/2 

f 



, {oiR - tti) (2 - /i) 
— +(-) 



2 (ai + aij)" r ((3 - m) /2) 
t +00. (22) 



Similarly to the case of fractional diffusion, we obtain the 
following particular time scale defined in dimensionless 
units: 



= max |l, ^4a^ai (l + + 2a^afl(l + AQ^ + 2a^afl(l + Afl^ + a^aj?))) + 2alal (b + 20n^ 

+2a^aR{5 + 20n'^ + a^afl(ll + 44f^2 + 4a^afl(3 + Un'^ + a^a,.)))) + (l + Ail'^) (2 + a^ai^ (4 + a^a,.(9 
+2af^ar (5 -I- 2a^aR 

)))) +2a^(l-f4r!2) a^(2 + 

aM"fl(4 + a^.aR^Q + 2a^aR (5 -I- 2a^aR) ))^ + (l + A^"^) 
xal(^^ + 2a^aR{Q + 2a^aR{lQ + a^aR{ll + Aa^aR)))^\ U'^aln'ala^iaL + aR) ) [>, (23) 
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where = T(i-'')/2/^r (2 - /i). 

We also notice that, for long collision time intervals, 
t ^ T, a reduction of the probability that two consecu- 
tive collisions occur. 



T2 < Ti, 



i)Tr 



< 



{t+T,r 



-00, 



is related to a decreasing of the parameter T. Thus, a de- 
creasing of the probability that two consecutive collisions 
occur at large time intervals, provokes a faster relaxation 
of the initially populated L levels. The inverse power law 
relaxations, occurring over long time scales, are acceler- 
ated by a decreasing of the probability that two consec- 
utive collisions occur at large time intervals. As in the 
previous case, this correspondence could be identified as 
an IZE over long time scales, occurring for large collision 
time intervals. 



C. Exponential memory kernel 

The case of exponential kernel, ^^{t) = Ae"'''*, is ar- 
gument of interest in stochastic processes [l^. The cor- 
responding distribution of collision times is 

sinh (tVl^ - 4, a) e"''*/2 
w^it) = 2A ^ , ^ . / , 7'>4A, (24) 



with finite mean time = •y/A. 

The asymptotic dynamics of the levels populations and 
the coherence term are obtained from the Laplace trans- 
form of the memory kernel, $7(1*) — A/ + u), and the 
corresponding behavior of Eq. (|B1[) in the limit w — 0^ . 
With this method, the following asymptotic forms are 
obtained: 



2nT^ {aL+ any r {-1/2) 



-3/2 



t. +00, 



(25) 



+ (-) 



("L - api) 



aL+aR 2{aL + any T (1/2) 



-1/2 



t — > +00, 



(26) 



describing the time evolutions of the coherence term and 
the whole L and R levels populations, respectively, over 
a long time scale, t^Tj, defined below. 



The convergence of the series expansion of Eq. (jBll) . 
leads to the following choice of the long time scale defined 
in dimensionless units: 




4a? q;3, (ttL 



rp3 

7 




8aia 



1 ^ 9 
aj^ 2aLOiR 

SaLttR 



9al 



j.5/2 

16a|Q;|j {uL + ur) 



rp3/2 

7 



(27) 



The analysis of the asymptotic dynamics reveals that an 
increase of the collision rate, i.e., a decrease of the mean 
time , provokes a faster relaxation of the initially pop- 
ulated level over long time scales, t ^ t^. This is a clear 
manifestation of the IZE over long time scales. 

For vanishing values of the mean time, T-y <^ 1, the 
long time scale tends to the value 
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ri°) = max < 1, — ^ 
1 (2^!)^ 



The IZE still reveals for t ^ T!f \ but the system ap- 
proaches the stable configuration (IT^ , this is compatible 
with the appearance of the QZE at earlier times. 



(28) 



D. Bi-exponential distribution 

We now consider the case of a bi-cxponential [2l| dis- 
tribution of collision times 

Wbeit) = PaDae-"-' + PbDbe-°'\ (29) 
with finite mean time 

PaDb + PbDa 



The — 



DaDh 



As in the previous cases, we plug the Laplace transform 
of the corresponding memory kernel 



DgDb + U (DgPg + j^fePfe) 

DaPb + DbPa + U 
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in the expression (jBl|) and we study the behavior for 
M — O"*". The Laplace inversion leads to the following 
asymptotic forms: 



t 

The 



t 

The 



2nTbe K+aflyr(-l/2) 

-3/2 

, t — > +CX), 



(30) 



-1/2 



{an - ul) 



{aL + anyv {1/2) 



t +00, 



(31) 



describing the dynamics over long a time scale, t » The, 
defined below. The convergence of the series expansion 
of Eq. ()B1| leads to the following particular time scale 
defined in dimensionless units: 



Ue = max < 1, ^ 



16a^/^ala% {aL + an) (^{a + bf + 4bn^{3a^ + 3ab + b^)^ 

+ 16 TtY'«'^'«iafl(« + b)^ (1 + 4f^') ((« + b) («! + a%) + SaaLafl) + 4 a^/^a^a^ K + an) 
X (a + 6)2 (1 + m^) (2 {aL + anf + 9aLafi) + ^iT^^^^ a^^^{a + b) {l + m^) ((a + b)ai + 5aalaR 
+ 10aia^(a + b) + 5aaLa% + {a + b)a%^ + 2 a^/^ j^2^(^^ ^ ^) (q,^ ^_ q,^) _^ 4f^2^ ^^,2 _^ ^2^^ 



X (2 (ai + aR) + tH^) 




2'^DlDln'ala% {aR + a^) 



r 



(32) 



where 



Dr,Dh 



+ Di,Pt, d^DaPb + DbPa 



The above analysis suggests that an increase of the 
collision rate, i.e., a decrease of the mean collision time 
Tfce, provokes a faster relaxation of the whole population 
of the L level, revealing the occurrence of the IZE over 
long time scales, The- 

In case of a vanishing mean time, T^e <C 1, the time 
scale for inverse power laws tends to the value 



(0) 

r^V = max. 



1 a[{a + bf + ibn^{3a^ + 3ab + b^)y 



(33) 



The IZE still reveals for t ^ T^'g\ but the system ap- 
proaches the stable configuration (fT2|) . This behavior is 
compatible with the appearance of the QZE at earlier 
times. 



The particular case Pa — 1 and Pb — gives the 
Poisson distribution of collision times, Wp{t) = e~*^'^° /tq, 
and it is described by the master equation ([S]), where 
Tfce = 1/-Da = To- Over long time scales, the dynamics is 
still described by the inverse power laws (PT|) . except for 
the different choice of the time scales, t ^ Tp, given by 
the following expression in dimensionless units: 



Tp = max-^ 1, -, 



1/2 







16Q;|a|j {aL + Ofi) 



+llaLaR {aL + aR) ) + 4to(9 (a^ 



SaLaR) 



- 4ro'/' {aL + aR) + 2tI 



(34) 
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For vanishing tq, the above time scale tends to the value 



r/^-* — max < 1, — , j 



(35) 



As in the previous cases, the IZE reveals for t » rp°\ but 
the system approaches the configuration of asymptotic 
equilibrium (1121) . Again, this behavior is compatible with 
the appearance of the QZE at earlier times. 



IV. CONCLUDING REMARKS 



time scales, are accelerated by a decrease of the proba- 
bility that two consecutive collisions occur at large time 
intervals. This correspondence could be identified as 
an IZE over long time scales, occurring for large colli- 
sion time intervals. The cases with finite mean collision 
times, show faster relaxation as the mean collision time 
decreases, revealing the occurrence of the IZE over esti- 
mated long time scales. 

For a vanishing mean collision time, the system ap- 
proaches the condition of asymptotic stable equilibrium, 
over the time scale the IZE reveals. This behavior is 
compatible with the emergence of QZE at earlier times. 



We consider a gas of identical chiral model molecules, 
interacting according to a Poisson distribution of collision 
times, described in Ref. 11]. We apply the celebrated 
prescriptions of the CTQRW to build up the master equa- 
tion driving the exact dynamics of the rotational internal 
energy levels of the model molecule for a general distri- 
bution of collision times. The resulting master equations 
turn out to be non-Markovian and the Lindblad form is 
recovered in the particular case of a Poisson distribution 
of collision times. 

The reasonable assumption that the molecule is far 
from resonant conditions and the dissociation energy, al- 
lows to consider as infinite the number of energy levels, 
this way, the asymptotic exact dynamics is analytically 
described. The whole L and R levels populations tend to 
the same asymptotic stable equilibrium configurations, 
for every random distribution of the collision time inter- 
vals. 

The levels populations and the decoherence terms ex- 
hibit inverse power law relaxations over estimated long 
time scales, for bi-exponential and power law distribu- 
tions of collision times, for fractional diffusion and in the 
case of exponential memory kernel. 

The cases with infinite mean collision times show that 
the inverse power law relaxations, occurring over long 



Appendix A: The General Master Equation in 
details 

This Appendix is dedicated to the study of the dy- 
namics described by Eq. ([SJ. Since the ground states 
are degenerate, the Hamiltonian H takes the following 
diagonal form: 

if = i?ijl_)(l_| 

+i?ijl+)(l+|+ 5] ^ £;„Jn,)(n.|, (Al) 

lU<,,)=7iti|l+)+7U|l-), (l+|l->=0, 

|1±) = Vfl) + si^^lU), i?i± = ± ^f^, 
(±) _ (±) _ (±) _ (±) _ 1 

this way, after long but straightforward algebra the time 
evolution of each element of the statistical density matrix 
is obtained. As follows, the frequency terms ar labeled as 
'-^v,r) — {El, — Eri) /H. We start from the evaluation of the 
master equation driving the dynamics of matrix element 



4 4 t 

i=LM.j=L,Rn=l m=l •' ^ 



where the non vanishing terms (/jJiitm/' (''') ^-re 
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Vfci'/' ir) = ^aL cos^ (Or) , ^(^^^f^) (r) = za^^ sin^ (Or) , ^^^-f^^ (r) = - (<^(^„-^/) (r)) * = - sin (2f2r) , 
Vfe'"^ (r) = aie-==-^-- cos (Or) , ^['^f^^ (r) = »a,.e-^-^-- cos (Or) , ^[['^^f;^ (r) = e^^-^-^ cos (Or) , 
41^.^' (r) = -f e--...^ = _2 (^^y ^ ^(U,2^.) = 3 ^,|sin(Or) e-^^.^ = 

^ {viti^Ur})* , ^(^^-';)(r) = «^sin(0r)e-3...-, ^^It^.^' (r) = -*aL, ^fc^/^ M = -^e<--^--^K, 
•P^it.'^ (^) = e-3.=.-ai, (r) = -aie— ^^^^^ cos (Or) . 



The term /Oii,,ii,(t) reads as follows: 
I 



i=L,R j=L,R n=l m=l ^ 



r 

and the non- vanishing terms ipnt^mp (r) are listed below: 



^t:uHT) = -alcos'{nr), ^[\%^^ir) = -alsinAnr), ^1^^^;^ (r) = -.^ sin (20r) , ^^^^^^r) = ai, 
y'bi'^ (^) = -^e-3.i.-cos(Or) = (<pi\^^i\^^ (r))* , ^(^^^^^^^ (r) = .a,. e^^.^.^ cos (Or) = (r))* , 

'^l'„';i"^T)=aisin(Or)e'-.^.-= (<^(i-^;)(r^^ (t) = sin (Or) e-3.i.- = 



r 



The term piii,ii{{t) is obtained by exchanging the sub- As regards the coherent term pi^iRit), we obtain 
scripts R and L. 



Pl.,lAt)=^^{Pl.ldt)-Plnln{t))+ Y.ilil r$(i-Ovi'tm;^(t-OPn.,m,(i')'^i', (A4) 

i=L,R j=L,R n=l m=l 

I 

with the following non-vanishing terms: / \_ (^/J + '^i) • 2 /o ^ 

(U,i«) ^ -»K+Q|t) sin(20r) ^(i;;^^ ) (^) = sin (Or) e^^^^fl 



_ ' (ua.) , . ^&3«"^ (r) = ^«?ee-3«..^ sin (Or) , 

n 1 ^ fy^ 'PuL (^) = «afle'"^«i«^ COS (Or) , 

^&if^(-) = -^H^cos(20r), W 
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^U,in^ (^) = -^e-3Hi«- COS (Or) , y'&i'/^ (r) = a^e^^^^^^ sin (fir) , v'fc^"^ W = aflaie<-=«^«-^^^^)^ 
^^i^d"^ {r) = -za^e--.^-- cos (Or) , v^fci^^^^ (r) = -|ai sin (Or) e— ^.i.-, 

^ilai"^ (^) = e— 3.1.- cos (fir) . 



The term P2l,2l{^) reads 



4 4 J 

P2.,2.(i)= Ullll I ^{t-t')vZ';i^Ht-t')pn,m,{t')dt', (A5) 

i=L,fl 7=L,iJ n=l m=l ^ 



with the following non vanishing terms: 
I 



'ft\u^ (r) = al cos (fir) , ^'-^^'^^^ (r) = t'^ sin (2fir) = (^f^^f^^^ (r)) * , ^f^^^'^ (r) = -^a^e'--^.- cos (fir) 
= {^t'l^^ {r))* , (r) = cos (fir) , (r) = ^ sin^ (fir) , <pf„-^;' (r) = -zai 

sin (fir) , ip^^^%'^ (r) = aie^-^^.i.r ^ ^(2.^2^.) ^ ^(2.,2.) ^ -aie— ^-^^^ sin (fir) , 

^t.i:\r) = -Sal <pt3:^{r)^^aLe^-^^■^^^, ^t'^^Hr) = -ale^'^^r..^^^ , ^i^^,^^ (r) = ale— 3.,i.- 
X cos (fir) , ^g'f^''> (r) = *aie— ^..i.- sin (fir) , <p^l%-^ (r) = -za^e— 3.,=.-, ^(2^.2-) (^) = _^e— ^-^-^ 



I 

The terms pi„,i„(t), Pih.i.(Oi P2R,2i<if)^ PsL^Shi^) sake of shortness, we write down only the term ps^^s^ (f), 

obtained by exchanging the indexes L and i?. For the 



P3.,3.(i)= E E EE r*(i-0v'S?4')(i-0Pn.m,(i'X (A6) 
i=L,R j=L,R n=l m=l ''^ 



and the corresponding non- vanishing terms, 



10 



(^) = ^ W) * = COS (Qr) , <pf^l^'^ (r) = e^^...- sin (17r) , ^^^^^^^^ (r) 

= <^K"^ (r) = al, ^1%-^ (r) = (<^f^-^^-' (r))* = ^^y^-) ^ ^-...,,3,.^ 

'/'^tf.'^ M = (^^i:^ M)* = ^t\t:^ (r) = -^^e— 3-^.^sin(f7r) , ^^^^^ (r) = -2ai, 



The terms pms,ms (Oi ^^^^ given by the following forms: 
I 



ms+2 ms+2 

Pm.,m.(i)= J2 J2 ^(t-t')vi7'"''Ht-t')pi,j{t')dt', 4<ms<N,-2, s = L,R, (A7) 



with the following non-vanishing terms: 
I 



2 

2 

- _^P»(Sm3+2-i=;n.,)r/ft (m^.m,) \ _ 2 j(i<;^3_i-E^,+i)r/ft (m^.m,) / \ _ _ iC^™, -S„.3 + i)r/ft 

2 

(m,,m,) _ „2 (ms,ms) / ^ _ _ ^ -Bm3+2)r/fi 



For the sake of shortness, we omit the master equations For times t ^ 1/^^,, the contribution of the oscillat- 
of the remaining terms, recovering the same structure as ing terms to the convolution is negligible and the above 
the above equalities. master equations get simplified forms: 
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x/ cO(t_t'){p,^ ,^(t')+p,^ ^Jt')_2p2,,2. (Olt^^', (AlO) 
Jo 

nt 

Pu,2h W ^ * (^2«i«Pu,2„ W - f^Pi„,2« W) + / - t') (i') + [t')} , (All) 



2 







2 /-t 

P2„2, (i) ~^J^<P{t- t') (i') + {t') - 4p2„2, (f) + 2p3.,3. (i')} rfi', 



(A12) 
(A13) 

Jo 

m = 3, . . . , Ns — 1, s — L, 

PN^,NAt)-(^l f '^{t-t'){pN^~i.N^,-i{t')- PN^,NAt')}dt', s^L,R. (A14) 

I 







The dynamics of the populations, pm^.m^ (0: ^^id the co- The reduced master equation provided by the adiabatic 

herent term, results to be decoupled from the time theorem [2^ reads 

evolution of remaining coherence terms, Pms,n fo'" 
every m,^,n.^i # 1 and s,s' = i,-R, undergoing "fast" 
oscillations |ll| and giving vanishing contribution for 

t > 1/rj. 



1. The adiabatic theorem 

We now prove the consistency of the above equation 
through the adiabatic theorem described in Ref. [28]. 



-p{t) ^-{[H, Pit)] + J\{t-t')!^-t [V, R{t,t')]-^{U [V'] ,R{t,t')}+U [VR {t, t')V]^ dt', (A15) 



with the following definitions and use full relations: 
I 



N, 



n[^] = A = QAQ+ ^ ^ P„^AP„^, Q = Pi^+Pi_, R{t,t')^eM'-'')[p{t')], Pi± = |1±)(1±| , 

s=L,Rn^=2 

Pn^ = \n,){ns\, ns = 2,...,N„ s^L,R, V ^ 0, U [V^] = aU P^^ + Pn, + 2 Pn.] , 

s=LM. \ n»=2 / 

alPi, + a'flPiH |1+>(1+I +Pi-,i- |l-)(l-l |1->(1+I |l+)(l-l , 
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For long timescales, i 1/fi, the oscillating terms give the following relations: 
a negligible contribution to the convolution and we get 



-\ {n [V^] ,R] = -^1^ {p^,,^, it') |1+)(1+| + (0 |1_)(1_|) + (f) + pi_a_ (0) 

x(|l+)(l_| + |l_)(l+|)- ^ al[2 J2 Pn..njn,){ns\+ PN^,N AN,) {N,\], 

s=L,R ns=2 

n = "'{/'2„2, (i') |1«)(1«| +pjv,-i,jv,-i {t') \Ns){Ns\ + +P3,A (0) |2.)(2.| 



s=L,R 



Tls— 3 



giving the following master equations 



(t') + Pi-4- (i')) +«i7i^^i"V2.,2, (0 



+«K7ij+^i V2h,2h (i') ) dt'. 



We are now equipped to build up the master equa- 
tions generated by Eq. (|A15|) . After some straightfor- 
ward algebra we recover the same equations describing 
the dynamics of the levels populations and off diagonal 
elements of the density matrix obtained in the previous 
Subsection. 



Appendix B: The populations of the ground levels 
and coherence 



This Appendix is devoted to a detailed study of the 
asymptotic dynamics of the populations and the coher- 
ence terms in the Laplace space, describing the behaviors 
for times t ^ 1/17. 

We consider Eqs. IM]), (EHJ, dSUl, 

and those obtained by exchanging the indexes L and i?, 
the corresponding Laplace transforms show an iterative 
structure of continued fraction difficult to inverted, thus, 
in order to evaluate pi^{t), pif.{t) and p'^{t), we assume 
that the number of internal energy levels is infinite. This 



way, Eq. (jA13[) gives the following difference equation 

Pn,+2{u) - 2 ( ^ + 1 ) Pns + l{u) + (") = 0, 



V2a2$(w) 
= 2, 3, . . . , s — L,R. 

whose solution is 

p„Au) = Mu) [x'A\u)y" + BAu) (x^^\u) 

Us — 2,3, ... , s — L,R 



where As{u) and Bs{u) are independent of rig and 



^ 2a2$(u) 
TT-s = 2, 3, . . . , s = L,R. 



±1 1 



2a2$(u) 



1, 



The Laplace transform of normalization constraint of 
the whole levels populations, J2s=l rJ^'^^ =iPn,{t) = 1, 
reads _r (") — Also, we consider 
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the case of functions positive for every u > 

and s = L,R, thus, A*^-* (u) > 1 for every u > 
and s = L,R. This way, the convergence of the series 

T,s=L,R T,Z=2 ^''(") (^+^ (")) giv'^s Asiu) = for ev- 
ery u > and s = L, i?, and we finaUy obtain 

P„^(u) = S,(u)(aL''(u))"\ = 2,3,..., s^L,R. 

The function Bs{u) is fixed by Eq. (jAl2p and the one 



obtained by the mutual exchange of L and R, 

for every s = L,i?. The dynamics of the populations of 
the lower energy levels, pi^ (m), pi^(u) and the coherence 
term, p'^{u), is describe by the matrix form M(u)-X(u) = 
Y, where M(m), X(u) and Y read 





1 


2 









a-. 



1 + 



[MH]3_3 = [M(u)]4_4 = 

X(^) 



(aL^^m) -2) -u^ 



«fl (!>(«)) 


2 


2 (a\Hu) 


(aL-)H 


-) 





[M(^.)]3 , [M(u)]4.3 = 







" 1 - 


PlL,li?.(") 


, Y = 






- Pi«aL(") - 




. 0. 



Finally, the Laplace transform of the coherence term is 
I 



while the population of the ground L reads 
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(u) = (u + 2ail>(w) + V^Fl(m)) |^2VI^ + 21]2) + Fr{u)) + al^{u) (sil^ + + u^/^Fr{u)^ j 
{u^ + m^) (V^ {y/n + FLiu)) (2u (V^ + FRiu)) + apH (5^^ + Fr{u))'^ 
+al^u) (V^(5V^ + Fi(M)) {V^ + Fr{u)) + 2al^u) {6V^ + Fl{u) + Fr{u))) ) J. (B2) 




The expression for pi^(w) is obtained from the above re- 
lation through the mutual exchange of the indexes L and 
R. We notice that the equality (|lip is fulfilled. The co- 
herence term described by Eq. (|B1[) . is evaluated for var- 
ious expressions of the memory kernel In each case, 
the inverse power law behavior of both p'^{t) and PL{t) 
emerges from the analysis of the corresponding Laplace 
transforms for m — > 0"*", while the long time scales are 
derived from the conditions for the convergence of the 
related series expansions. 
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